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Abstract 

The stopping power of a metal for a slow structureless impurity varies linearly with projectile 
velocity. We show that the coefficient of this linear behavior (friction coefficient) is determined ex- 
actly in a static ensemble Kohn-Sham scheme, by extension of the Shifted Fermi Surface procedure 
originally derived for an homogeneous jellium. We prove that the friction coefficient is determined 
only by local characteristics of the system. The error incurred when adding a spurious non-local 
contribution is illustrated for a simple ID model. 
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I. INTRODUCTION 



The stopping power of a free (homogeneous) electron gas (FEG) for a structureless atom 
or ion varies linearly with projectile velocity in the limit of low velocities (see, e.g., Ref. 
This result has been proved at the level of both linear and non-linear theories of stopping. 
The coefficient of this linear behavior is referred to as the friction coefficient and the asso- 
ciated process as electronic friction. Such a behavior has been observed in experiments for 
ion stopping in metals that are expected to mimic a FEG in this respect. Examples are for 
H, D and He ions in Al^, H and He ions in Au and Cu^ or H and He ions in Al, Zn and Au^. 
It has also been observed in grazing collisions of ions on surfaces, like H and He ions on the 
Al(lll) surface^. 

One important field in which electronic friction might play an important role, is that of 
surface chemistry, as a candidate for energy dissipation by adsorbates. It is crucial to assess 
its importance for such processes as atomic and molecular adsorption, diffusion on surfaces, 
range of hot atoms following molecular dissociation at surfaces, etc. Then, the adsorbate is 
immersed into a strongly non homogeneous electronic density which raises questions as to 
the relevance of values of friction coefficients obtained from calculations based on the FEG 
model. It is, therefore, desirable to develop methodologies that can make use of present day 
band structure codes to evaluate the friction coefficient. Such a task has been endeavored by 
Trail et alP, based on the formulation proposed by Hellsing and Persson~. Results have been 
obtained for H and D on Cu(lll)^^ or for H2 on Cu(lll) and N2 on Ru(OOOl)^. However, 
some questions are raised by the latter results. Firstly, the H/Cu results are divergent at a 
certain value of the atom-surface distance. Secondly, the friction coefficient for N2/RU may 
reach values an order of magnitude larger than usually obtained from the FEG model. 

In the present work, we reexamine the formulation of the friction coefficient determination 
from first principles in order to ensure that we get a consistent first order approximation 
in the projectile velocity. Our approach is based on Density Functional Theory (DFT) and 
extends the methodology developed for the homogeneous FEG to the inhomogeneous case. 

Atomic units are used throughout. 
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II. ADIABATIC THEORY 



We consider a medium composed of electrons moving in a periodic array of fixed nuclei. 
We are interested in the projectile energy loss due to the electronic state perturbation by the 
projectile motion, the role of the projectile frozen-lattice interaction (i.e., with the lattice 
nuclei) being trivial. The "external" potential, in which the electrons move, is composed of 
two parts: one, l^iat, due to the fixed nuclei and the other, Vepit), to the moving, structure- 
less, projectile. As the projectile position with respect to the lattice, Rp, is time-dependent, 
the friction coefficient is a function of Rp. We suppose that the static problem, with the 
projectile at rest, has already been solved and restrict our study to the limit when the pro- 
jectile velocity goes to zero, since the friction coefficient is associated with the first order 
approximation to the stopping in its power series expansion as a function of the projectile 
velocity. In principle, the density of the system, as a function of time, may be determined 
through the time- dependent Kohn-Sham (KS) scheme of the Time Dependent Density Func- 
tional Theory— (TDDFT, for a general presentation, see e.g., Ref. llSl ). However, in the low 
velocity limit, it is not necessary to resort to TDDFT as we show below. 

Our strategy is as follows. In the present section we study the behavior of the system in 
the adiabatic limit. In the next section we show that the non-adiabatic corrections do not 
contribute to the friction coefficient. 

Let us start with the projectile at rest in the lattice frame without interaction with 
the electrons. We could first solve the static problem by turning on the projectile-electron 
interaction and, subsequently, set the slow projectile motion. However, in the adiabatic 
limit, we may as well first set the projectile motion, without interaction with the electrons, 
and determine the transformation from the lattice to the projectile frame. We then turn 
on the projectile-electron interaction adiabatically. The latter procedure is justified since in 
the adiabatic limit the final state is independent of the evolution path. 



A. Transformation to the projectile frame 

We first consider the projectile moving without interaction with the electrons. We assume 
that the projectile velocity, v, is constant so that the transformation from the lattice to the 
projectile frame is Galilean. Let be the position vector of electron i with respect to an 
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origin fixed in the lattice frame and rf w.r.t. tlie projectile. The N-electron wave-function 
\E'T(ri, rAf) becomes in the projectile frame (discarding in this section, for short, the 
energy phases that play here no role): 

vl/p(rf , Rp{t)) = e-'Sf-^-< ^T{r„ r^) (1) 

with Vi = rf + Rp{t). The density in the target frame, and in the absence of projectile- 
target interaction, is exactly known through a static Kohn-Sham scheme. This static Kohn- 
Sham scheme is transformed into the projectile frame by a similar Galilean transformation 
involving the multiplication of each KS orbital in the lattice frame, iff , by exp{— i?; ■ r^}. 
As the KS orbitals are Bloch functions: 

<,(r^; Rp{t)) = e—^'^^Ur) = e^^-e^^^-^)- «.,(r) (2) 

where Uk,i is a periodic functional. In the absence of a gap close to the Fermi surface: 

= ^l-Ur) + Oiv) (3) 

When substituting v^^_^,j for ip^^ (disregarding the common time dependent phase factor 
exp{ii> -iip}), we introduce a shifted KS scheme in the target frame, i.e. a Fermi distribution 
shifted by ~v (SFS). It can be shown that this shifted KS scheme is correct to first order in 
the velocity. Firstly, the total number of states within the shifted Fermi surface is conserved 
up to first order in v. Indeed, for vanishingly small v, and in the absence of a gap in the 
vicinity of the Fermi surface: 

j dk p{k -v) = J dk p{k) + j dkpi-v ■ kp) pikp) + 0(^7^) (4) 

where p is the density of levels and the volume integrals are over the volume bounded by 
the unshifted Fermi surface. Since the Fermi surface is symmetric with respect to the origin 
of the Brillouin zone, the quantity {—v ■ kp) dkp is exactly compensated by the same term 
for k = —kp. Secondly, the density is also exact to first order in v. 

n^'ir) = [ dk\^l_,^,{r) + e'^''-''^-^v-VkUk,{r)\' + 0{v') 
Jfs 

= [ dk\\^l^,,,ir)\' + 2^[uUr)v-V^u^4r)]]+0iv') (5) 
Jfs ^ ^ 
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The symmetry with respect to the origin of the Brillouin zone imphes that ^ = U-k,i so 
that the second term in the integral gives a zero contribution when integrated over k, which 
yields: 

n{r)= I dk \vUr)\' + 0{v') (6) 

In other terms, the SFS prescription gives the correct density up to first order in v which 
means that the corresponding KS scheme is exact to first order in the velocity. We may 
re- write the previous expression for short, as: 

n^(r^; Rp{t)) = n^,(r^; Rp{t)) = ^ \^J{r)\' + 0{v') (7) 

i 

where Fi corresponds to the Fermi distribution shifted by —v. 

As expression ([7]) shows, the time dependent KS solution in the projectile frame can 
be formulated in terms of an ensemble KS procedure in the lattice frame. When v kp 
and for a given orientation of fc, the occupation of orbital yjfc is 1 for k < and zero for 
k > km, where km = kpk — v. The Fermi surface being a surface of constant energy, this 
entails that the population decreases as the energy increases as required for the validity of 
the Hohenberg-Kohn theorem for ensembles^^. Then, the latter theorem proves, in the non- 
degenerate case, that the exact state of the system in the projectile frame is also described 
by an ensemble in the lattice frame, univocally defined from the K-S ensemble since the 
latter determines both the ensemble and the density. 

The system we are considering here involves degeneracies, so its ensemble is not deter- 
mined univocally by the density. However, the energy remains a functional of the density 
and any quantity that can be expressed in terms of the density is also univocally defined. 
These properties are the only ones required for the validity of the following discussion. 

As a conclusion, the Galilean transformation from the target to the projectile frame of the 
exact N-electron state yields, to first order in the velocity, an ensemble of states which are 
defined in the (static) target frame. In other terms, we have transformed the time-dependent 
problem in the projectile frame into a static one in the target frame, at the expense of using 
an ensemble defined by the Shifted Fermi Surface (SFS). 
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B. Projectile-target interaction in the adiabatic limit 



1. Adiabatic evolution 

Let "^jiti) be the exact solution of a time dependent problem when the system is in the 
eigenstate j of the Hamiltonian at the initial time t = U. 

H{U)^j{U) = Ej{U)^>,{U) (8) 

If the evolution is adiabatic, the system remains in the eigenstate j of H{t) at any time, 
i.e., it evolves in such a way that the wave function is at any time the solution of a static 
problem. Evolution from two different initial states cannot lead to the same state at time 
t: adiabatic evolution is reversible. For simplicity, the derivation below assumes that the 
eigenenergies are not degenerate--. Accordingly, the function j{t) is determined at any 
time t (save for a time-dependent phase factor) by the energy Ej(t), i.e., the evolution of the 
system state is determined by the knowledge of dEj (t) / dt. This one to one correspondence 
between energy and state is central to our derivation below. 

The same correspondence prevails for the evolution of an ensemble, which may be defined 
by the density operator: 

j 

where the weights pj are constant in time. The energy is given by: 

Eit) = Trim Hit)} = J2 P3 Ejit) (10) 

j 

Again, the system is at all times characterized by quantities that are the solution of a static 
problem. Furthermore, when the p/s satisfy the same condition as that for the validity of 
the HK theorem for ensemble, the ensemble is determined at any time by its energy. 

2. Application to the present problem 

In the present problem, the adiabatic evolution is associated with that of the projectile- 
lattice relative position, i.e. with the variable Rp. What remains to be derived is the KS 
scheme associated with this adiabatic evolution. To this end, we switch on adiabatically the 
projectile-target interaction by writing the external potential as: 

VeAr, t) = X{t) Vep[r, Rp{t)] + V,,t (H) 
6 



The function X{t) is an arbitrary function that varies from to 1 in the time interval [ti,tf]. 
Saying otherwise, we introduce through the function A(t) a fictitious evolution starting from 
an initial state where the projectile and electrons do not interact. Although we could keep 
Rp fixed in this process, we conserve in flTT]) its time- dependence. As can be seen below, the 
conclusions of the present section are independent of the presence or absence of the latter 
time-dependence. We assume that it is always possible to choose the arbitrary function 
A(t) so that dX/dt is small enough, at any time, to ensure adiabaticity. The adiabatic 
approximation corresponds to a zero order approximation in dX/dt (see,e.g.,— ). 

Our aim is to determine the KS orbitals, at each time, when the system evolves adiabat- 
ically. To this end, we calculate the energy variation in the interval dt: 

drn{r,t) — V^p + X{t) v V R^Vep\ (12) 

where n is the exact density at t. According to the discussion of section Hi B 11 the KS orbitals 
must be such that they give the exact energy variation ( IT2l) in the limit of vanishing dX/dt. 
We suppose that, at a given time we have determined the KS orbitals for the ensemble 
defined by the distribution Gi and study the time evolution in the interval + dt]. 

Consider the orbitals (p'^{r,t) evolving adiabatically from the KS orbitals at t° for the same 
ensemble, i.e., assuming that Gi is independent of time: 

n'^{r,t) = Y,G.\^1ir,t)f (13) 

i 

where the orbitals are the solution of a static problem 

Hj,s^tir,t) = st{t)iptir,t) (14) 

with 

Hks = + X{t) V.pir; Rp) + t;„(r; K(r, t)]) (15) 

and Vn includes the Hartree potential, the lattice potential and the exchange correlation 
potential of the static problem for the KS ensemble defined by the occupations Gi. The 
question is now: what is the relation between these adiabatic KS orbitals and the KS orbitals 
associated with the exact adiabatic evolution of the system? Introduce the energy functional: 

Ev[n] = ^Gi£i{t;[n]) - j drn{r,t)vy,c{r;[n{r,t)]) 




, , , n(r,t)n(r',t) ^ . , 



(16) 



The functional Ey[n\ gives the exact energy of the static ensemble with density n, when the 
ej's, the exchange-correlation potential Vxc and the energy Exc are those of the static KS 
ensemble defined by the occupations Gj. At t^, n = n"" = and -Et,[r;.°] = E{t^) where 
and E{t^) are the exact density and energy. The question is: what is the relation between 
i?^,[?7,"] and E{t) at t = t° + dt? We split the variation of over the interval dt into two 

terms. The first one corresponds to the variation of ra": 



dt 



dr 



d£;Jn1 



dn" 



iO(r,tO) 



dn''(r,t) 
dt 



(17) 



This term is identically zero because at t = t° the energy functional is stationary around the 
exact density n^{r,t^). The second term corresponds to the explicit dependence on A(t) V(.p 
in f|T6|) for a fixed n"". 



AE^ln"] 



dt 



deyt) 
dt 



j drn°(r,t°) 



dA 
dt 



Kp + A(t) V ■ VhpKp 



(18) 



which is identical to ( !T2|) . As a conclusion, at t" + (it the y?" orbitals are still the KS orbitals 
for the exact adiabatic evolution when the ensemble is defined by the same time independent 
distribution Gj. This result may be applied to the evolution over the whole interval [ti,t/]. 
Now, for t = tj, i.e., in the absence of projectile target interaction (A(t) = 0), the distribution 
Gi is exactly given by the SFS prescription, Gi = Fi,'^i. This means that the solution of 
the problem at t = tj (A(t) = 1) is also given by the SFS prescription—. 

We are led to the following conclusion: the density, in the adiabatic limit, can be deter- 
mined, at each position of the projectile, using a static ensemble KS scheme in which both 
the projectile and lattice potential are fixed. The fact that the energy is stationary with 
respect to an arbitrary variation of the density around its exact value (see Eq. [T7|) plays a 
central role in reaching this conclusion. The SFS Kohn-Sham scheme is an exact scheme for 
the determination of the system density for an impurity traveling through an inhomogeneous 
medium under adiabatic conditions. 
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III. FRICTION COEFFICIENT 



The force acting on the projectile is: 

Vr^E{Rp) = - j drn{r- Rp) V RpVep{r- Rp) (19) 

which involves only the derivative of Vep since, as for (fT7|) . the derivative of the energy 
functional with respect to density is zero. 

When the density in (IT^ is the static density no(r; Rp), the force is the static force, i.e. 
the limit when the projectile velocity v goes to zero. When the projectile moves with a small 
but finite velocity t?, the associated density is nv{r\ Rp) and a dissipative process takes place 
which corresponds to friction. The friction coefficient for motion along the direction v is 
defined by: 

J-* = lim - j / dr/\n{r- Rp) v ■ V R^VeP^r; Rp) \ (20) 

where 

An{r-Rp) = n,{r- Rp) - no{r; Rp) (21) 



A. Non-adiabatic corrections 



When is calculated in the adiabatic SFS approximation of section Hi Bt the dependence 
of n„ on V arises entirely from the SFS. However, we have to check whether non-adiabatic 
corrections can be neglected when is calculated to first order in v. 

Deviations from the adiabatic approximation correspond to inelastic transitions between 
adiabatic states due to the variation of Rp{t). Let us evaluate the contribution of these 
inelastic transitions to (fT9l) when Rp{t) varies from R^ to Ri. We study the evolution of a 
KS orbital equal to the adiabatic orbital v^fco.i,, for Rp(to) = Rq. We set v ■ Rp(t) = Z = vt 
and, to simplify our notations, hereafter only mention explicitly the dependence on Z. For 
Z > Zq, the KS orbital becomes i'koMy which we express as: 

ipko,ioi^) = y~l'^'^''j(^) Vk',j{Z) exp I - - / dZ'ek',j{Z')\ (22) 

with the initial condition ipko^ioiZo) = ipko,io{ZQ). Using the time- dependent Schrodinger 
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equation, we obtain^^: 

^Cfc/j(Z) = -^Ck,iiZ){(pk',jiZ)\-^ipk,iiZ)) exp - ^ I dZ'[ek,i{Z') - e^'jiZ')] | 

k,i JZo 

(23) 

A first order approximation yields: 

Cfc/j(Zi) = 5k'jMio-J dZ{ipk',jiZ)\^ipk^^ioiZ)) exp j dZ' [eko,ioiZ') - Ek'jiZ')] | 

(24) 

Now, the V dependence of tlie excitations (i.e., tlie non adiabatic contributions) is entirely 
governed by the exponential terms in f l23|) or fl24|) . They oscillate rapidly for vanishing v, 
which quenches dramatically the transition probability when the two states are non degen- 
erate. Then, as is well known, the transition probability does not vary as a power low in 
V and is certainly not linear in v. This leads to a vanishing contribution to the Friction 
coefficient. However this is not the case for degenerate states or near-degenerate states: the 
transition probability is not quenched by the oscillations caused by the exponential if the 
energy difference between two states is of the order of v. We have to face this situation in 
our problem, since the occupied states belong to continua. We must, therefore, evaluate the 
corresponding contribution to the force (fT9l) at Zi. 

dVep{r; Z) 



S[Z,) = ^E{Z) 



Z=Zi 



J dr [nnonad(^'; ^l) - ^ad(^'; ^l)] 



dZ 



(25) 

z=z. 



where 



^nonad(r"; ^l) = / dfc (^l) |^ 

Jfv 

nad(r;Zi) = / dk\ifk,^iZ,)\^ (26) 
Jfv 

and the integration is inside the unshifted Fermi surface (FV). Using ( l24l) and keeping only 
the first order term: 

dVep{r; Z) 



S{Z,) = IZy^^dfc j dfc'23?[(<^fc,(Zi)| 



dZ 



l^k'jiZi)) 



Z=Zi 

j dZ{ipk'^j{Z)\-^ipk,i{Z)) exp I - -J- y dZ' [ek,i - Sk'j] } (27) 



Remembering that ipk.,i is an eigenfunction of the static {Z fixed) Schrodinger equation, we 
obtain readily: 

{Vk,i{Z)\^^^^^^^\Vk',j{Z)) = {ek,i - ey,j){^kAZ)\-^Vk',j{Z)) (28) 
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The matrix element in the r.h.s. depends only on properties of the adiabatic functions ipk,i 
and not on the velocity. However, from the preceding discussion, we know that inelastic 
contributions are only appreciable if {sk^i — Ek'j) is of the order of v. Furthermore, the 
integration over k' in fl27|) for a given k can be transformed into an integration over e^'j- 
The transformation from one integration variable to the other involves only properties of 
the static electronic structure and is independent of the velocity. Again, the range of the 
integration over Sk'j is of order v around ek,i = £k',j- As a consequence, the quantity S is 
of order v"^. 

Until now we have only considered the first order approximation to the Z-dependent 
problem as defined in (^^. However, it can be easily verified that each increase in the 
perturbative order introduces an additional integration over the energy of intermediate states 
and, therefore, an additional factor of v. As a consequence, the non-adiabatic contributions 
to the force on the projectile are of order v"^ at least, which means that they do not contribute 
to the friction coefficient. 

We conclude that both no and (to first order in v) can be determined exactly by the 
adiabatic SFS-KS scheme. So the latter scheme provides an exact procedure to determine 
the friction coefficient. 



B. Corollary 



The SFS-KS scheme relies only on local properties of the system, i.e., the determination 
of the density only requires information on the electronic state for a fixed value of Rp. 
It does not require information on the variation of any electronic quantity with Rp. Conse- 
quently, any alternative to fl20|) for the evaluation of the friction coefficient must satisfy the 
same condition. This provides a powerful tool to evaluate the validity of procedures aiming 
at the evaluation of Ti,. 



We may apply this condition to analyze the procedure of Trail et al^^^ (see also Ref.ll3l). 
The latter authors use the expression of the friction coefficient proposed by Hellsing and 



Perssoni^ (see also Ref. 
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)• 



Jkp = 2nkl dkp / dk' 



F -F 
,2 

V ■ Vrp\^ VKs{r,Rp] [n]) 



(29) 
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where vks is the full KS potential, kp the Fermi momentum and ip^ (resp. ip~) satisfies 
outgoing (resp. ingoing) boundary conditions. From the derivation of fl29|) in Ref. |lO| (or the 
alternative derivation in Ref. Il2l). it is not clear that it consists in a first order approximation 
in V, though when the impurity is moving through an homogeneous jellium, it can be proved^ 
that is equivalent to ( 12U]) to first order in t>. A key property, in the latter proof, is the 
invariance by translation of the free electron gas state in the absence of an external potential. 
For the inhomogeneous case such an equivalence has not been established. At any rate, in 
the latter case, and in view of the previous discussion, it is clear that the derivative of vks 
with respect to Rp in fl2^ must be carried out for a constant density, otherwise it would 
introduce a non-local contribution (associated with the variation of n when the projectile 



moves). However, the authors of Ref. Iq 
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- 13 use a finite difference method to calculate the 



derivative of vks'- they determine the KS potential for the projectile at rest and for two 
different values of -Rp, say -Rp + 6Rp and Rp — 6Rp. In so doing, they include a term 
associated with the variation of the function v^s with -Rp. The latter term, being non-local, 
introduces an error in the evaluation of the friction coefficient. 



IV. ILLUSTRATION WITH A ID MODEL 

The previous discussion has raised two questions. Firstly, how important is the error due 
to the non-local term in the evaluation of the friction coefficient? Secondly, is expression 
(!29|) valid for the inhomogeneous case? To answer these questions, we have built a simple 
ID model amenable to an exact numerical solution. Full details are given in the Ancillary 
materiali^. As we are concerned with general principles, our model needs not represent a real 
system, although we keep speaking of "electrons". Though units are here irrelevant, since 
we are only interested in relative results, we keep using atomic units so that the "electrons" 
have a mass of 1. The interaction between two particles located at z and z' is described by 
a screened regularized coulomb potential: 

vf{z,z') = ±e-'^\^-'y[\z-z'\+e] (30) 

{v^ between electrons and between electrons and "positive" particles). In the numerical 
application, we have used a = 2 and e = 0.1. We solve the problem at the level of the 
Hartree approximation. The Fermi momentum for a ID paramagnetic homogeneous jellium 
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FIG. 1: Sketch of the background potential used to create an inhomogeneous jehium. 

with linear density n is = n7r/2. We introduce the Wigner-Seitz distance Zs = 1/n, i.e., 
the length of the interval enclosing one charge. We create an inhomogeneous jellium by 
introducing a background step potential defined by: 

^'step(^) = 0.5AV{1 + COs[7t{z-Zi)/{z2-Zi)]} {Zi < Z < Z2) 

= AV {z< Zl) 

= {z> Z2) (31) 

We use throughout Zi = —0.5 and Z2 = 0.5. The value of AV {AV < 0) is fixed by imposing 
values of zl and z^ at z ^ — oo and z ^ oo respectively. The Fermi energy, ep, is then 
given by = (fcf )^/2 + = (fc|')^/2, so that rii > n2 and zl < z^. The energy diagram 
is represented in Fig. [H We introduce a background of positive particles such that the 
Hartree potential is zero in the absence of impurity (i.e., the density of positive particles is 
everywhere equal to that of the electrons). As an example, we plot in Fig. |2]the density 
for zl = 1, Zg = 2 and zl = l,Zg = 1000. In the latter case, the behavior of the density is 
qualitatively similar to that of a surface. 

The interaction between an "electron" at z and the impurity at Zp is represented by 
the potential v~{z, Zp) with a = 4 and e = 0.1. This potential supports one bound singly 
occupied state with energy -1.275. For all results given below, the projectile is located at 
Zp = 0. Calculation of the density in the presence of the impurity (even for v = 0) requires 
to confine the system into a box because the Hartree potential behaves as sin(2A;f2^ + 7i,2)/^ 
far from the impurity. We have used a finite interval z G [—5, 5]. This does not invalidate our 
conclusions since we are interested in the difference between two expressions of the friction 
coefficient for a given system. We have checked that our conclusions are independent of the 
box size and the implementation of the contour conditions. 

13 



0.8 ■ 



S 0.6 



a; 



0.4 



0.2 ■ 



z = 1 

S 

z^ = 2 



-30 -20 



-10 10 
Distance 



(a) - 



20 30 



— 05 
OS 



1 . 1 . 

a/wwwwaaaAi 

11 

- 


1 , 1 
(b) 

- 


z^ = l 




S 




X = 1000 




s 







-30 -20 



-10 10 
Distance 



20 30 



FIG. 2: Density in the presence of the background potential (j3ip with zi = —0.5, Z2 = 0.5, zl = 1 
and (a) Zg = 2 or (b) z^ = 1000. The density has been divided by the density of a uniform jehium 
with Zg = 1. 

In a first step, we assume that the impurity, the step potential and the associated back- 
ground of positive charges move together against a uniform jellium with Zs = zl- The friction 
coefficient can then be calculated in three different ways: 



(i) From the ID expression equivalent to fl2U]) : 



(0 



lim - 



dz [n^{z) - riQ^z)] 



Vep{z - Zp) 

+ v^t^^{z - Zp) + v+{z- Zp) 



(32) 



where Ustep and u+ are respectively the step potential and potential due to the back- 
ground of positive charges. 

(ii) It is easily shown^^ that it may be also calculated from the energy loss or gain associ- 
ated with reflection and transmission by the potential. When v < kJ^'- 



(33) 



where is the reflexion probability at the Fermi level for the static case. This very 
simple expression gives us a good check on the accuracy of our SFS calculations. 
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TABLE I: Friction coefficient when the impurity moves together with the background potential. 
Results in the first two columns are for double occupancy of the bound state and in the last two 
columns for single occupancy. 
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0.75 ; 1 


1 ; 1.5 


1 ; 2 


1 ; 5 


(Eq.Ell) 


0.576 


0.472 


0.545 


1.06 


(Eq.El 


0.573 


0.466 


0.546 


1.07 




0.576 


0.476 


0.560 


1.12 



(ill) From the ID version of (l29|l : 

= 2n J dk j dk' J dz[(pk'iz)]* ipk{z) 

^ VKsiz- Zp;[n{z,Zp)]) 5{eF - Sk) S^Sp - Sk') (34) 



dZp 

As explained above, the derivative with respect to Zp is evaluated for a constant 
density n, i.e., the quantities entering flM|) must be local. 



The equivalence between fl51|) and fl5^ or ( 1551) is due to the fact that, under the present 
conditions, the transformation of the system state from one position of the full external 



potential to another one involves merely a translation (see Ref. Il9| for more details). 

In Table [H we summarize the results obtained for some values of (z];z^). Results are 
the same for the three expressions of the friction coefficient within the accuracy of the 
calculations. This gives us confidence in our numerical procedures. 

We now consider the case of real interest for our model: that of the impurity moving 
alone, against the step-potential and background of positive charges, i.e., an impurity moving 
within an inhomogeneous jellium. The calculations only differ from the previous ones by 
dropping t>step and f+ in (!32l) and noting that the dependence of vks on Zp in for a 
constant density, arises entirely from Vep. Also, f p3|) is no longer valid. Results are given in 
Table M 

The first observation is that fl5^ is no longer equivalent to fl52]) . Therefore, the expression 
of Hellsing and Persson does not provide the correct first order in v for the energy loss of 
an impurity in an inhomogeneous medium. The error may be quite appreciable, as shown 
by our model. Note that the trivial difference mentioned above between the expressions 
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TABLE II: Same as Table |T] when the impurity is moving with respect to the fixed background 
potential. The last line gives the results obtained by inclusion in ()34p of the spurious non-local 
term when the derivative of the Kohn-Sham potential is calculated along (|35p . 





0.75 ; 1 


1 ; 1.5 


1 ; 2 


1 ; 5 


From Eq. [32] 


0.458 


0.321 


0.440 


0.71 


From Eq. [Ml 


1.944 


2.36 


2.59 


4.41 


Ad from jMl) in ^ 


0.404 


0.254 


0.379 


0.484 



for a moving and fixed background potential precludes a numerical error as the cause of 
this discrepancy. The fact that our ID model may be qualified as unrealistic cannot infirm 
the relevance of our conclusion since a single counterexample is sufficient to disprove an 
assertion. 

We evaluate now the error incurred when taking into account the change in the potential 
function in the expression of the potential derivative, as done in the finite difference proce- 
dure of Trail et al^->^^^. In place of the derivative of the KS potential for a constant density 
in f lMj) we use the expression: 

Av = {vKsiz, Zp + h; [nzp+h]) - Vks{z, Zp - h; [nzp^h])} /2/i (35) 

In actual calculations we have used h = 0.01. The two KS potentials in f l35|) are the result 
of a calculation where the projectile is at Zp + h and Zp — h respectively. For that reason, 
( !35l) includes a non-local contribution coming from the modification of the potential function 
when moving the projectile from Zp + h to Zp — h. Results are given on the third line of 
Table [Til They differ strongly from the previous ones (second line in Table [IT|) . up to nearly 
an order of magnitude. This difference is not due to the simple algorithm used in calculating 
the derivative as in ( 135|1 . If we use the same algorithm in expression ( IMj) . while keeping the 
density constant, results agree with those given on the second line of Table [ITl Changing h 
(within reasonable bounds) does not change either the conclusion. This demonstrates the 
significant error introduced by the functional derivative of the Kohn-Sham potential with 
respect to density. Note that the third line of Table [Tll should not be compared with the 
first one since it corresponds to an incorrect calculation of expression ( IMj) . the latter being 
intrinsically incorrect for an inhomogeneous system, as verified above. The two errors being 
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totally unrelated, there is no reason why they should compensate each other. The fact that 
the results of the third line in Table HT] are closer to those of the first one must be considered 
as merely accidental. 

V. CONCLUSION 

Starting from first principles, we have shown that the friction coefficient is, in principle, 
exactly determined by an ensemble Kohn-Sham procedure, the ensemble being defined by 
a Shifted Fermi Surface. We have shown that it depends only on local properties of the 
impurity /tar get system in the sense that its determination does not require information on 
the variation of the system density with impurity position. It may be of interest to check 
whether this conclusion could also be reached for other applications in which the low energy 
behavior is determined through an adiabatic approach. 

To our knowledge, this important constraint has been overlooked. As a consequence 
erroneous values of the friction coefficient have been obtained. This is the case, for example, 
in the work of Trail et al.^. In addition, the latter authors use expression (!29|) of Hellsing and 
Persson for the friction coefficient, which, though correct for an impurity moving through 
an homogeneous system, is not valid for the inhomogeneous case. 

That the error may be dramatic is illustrated by the divergence found by Trail et al.—^. 
They show, in their analysis, that the divergence is caused by the derivative of the system 
spin with respect to Rp. The latter quantity is basically non-local (locality being used 
here with the meaning defined above) and, therefore, the associated contribution to friction 
is spurious. The latter authors wrongly attribute the divergence to a breakdown of the 
adiabatic approximation. However, if the evolution for f — )■ is not adiabatic, this means 
that there is a discontinuity in the evolution of the system. In the case studied by Trail et 
al, no such discontinuity exists. So, the divergence can only be the result of an error in the 
evaluation of the friction coefficient, as confirmed by our analysis. 

Another conclusion can be drawn with respect to the evaluation by Luntz et ai— of 
the "local approximation for friction" (LDAF)^. In the latter approximation, the friction 
coefficient is evaluated as a weighted average of the friction coefficient for the impurity in 
a homogeneous electron gas having the local density at each point of the inhomogeneous 
system. The calculations of Luntz et alM, based on and are incorrect for the 
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inhomogeneous case and, therefore, the difference they find with respect to the LDAF cannot 
be considered as an evaluation of the error incurred when using the LDAF. 

Finally, we may remark that the evaluation of friction with present day band structure 
codes requires an SFS calculation since, up to now, no well founded alternative exists. 
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